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Abstract 

We present new analytical self-similar solutions describing a collapse of a 
massless scalar field in scalar-tensor theories. The solutions exhibit a type of 
critical behavior. The black hole mass for the near critical evolution is ana- 
lytically obtained for several scalar-tensor theories and the critical exponent is 
calculated. Within the framework of the analytical models we consider it is 
found that the black hole mass law for some scalar-tensor theories is of the form 
Mbh = f(p~Pcr)(p—Pcr)' y which is slightly different from the general relativistic 
law Mbh = const(p — p cr ) 7 . 

1 Introduction 

Critical phenomena in gravitational collapse have attracted much interest since the 
pioneering work by Choptuik [I] and they are considered to be one of the greatest 
successes of the numerical relativity. The essence of the critical phenomena consists 
in the fact that, just at the threshold of the black hole formation, the field dynamics 
becomes simple and exhibits discrete or continuous self-similarity despite the nature 
and non-linearity of the collapsing matter. The critical solution separates the solutions 
with a black hole formed from those without a formation of a black hole. The mass of 
a black hole formed in near critical collapse obeys the power law 



M BH = c( P - Pcr y (i) 

where the parameter p describes the strength of the initial data and p cr is its 
critical value. The exponent 7 is, in general, matter dependent and C is a constant 
which depends on the initial field configuration. For detail discussion of the critical 
phenomena we refer the reader to the Gundlach's review [2j and references therein. 

Parallel to the numerical investigations , an analytical approach has been under- 
taken in order to understand more deeply the critical collapse. Although fascinating, 
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this task was proved to be very difficult. The analytical approach is based on ap- 
proximate (perturbative) techniques or exact toy models which cannot be considered 
as completely realistic and satisfactory. Nevertheless, the analytical methods, and es- 
pecially the exact solutions provide valuable independent insight into the subject. It 
was realized by Brady [3] and Oshiro, Nakamura and Tomimatsu |5] that the one pa- 
rameter family of exact self-similar real massless scalar-field solutions first discovered 
by Roberts 1 |3| can be applied to explain some basic features of the critical collapse. 
Solutions from this family with p > 1 describe the formation of a black hole and to 
the leading order around the critical value p cr = 1 the black hole mass turns out to be 
M ~ (p — l) 1 / 2 . The critical exponent in this toy model is therefore 7 = 1/2 and is 
slightly different from the numerically calculated value 7 = 0.37. In order to gain as 
much as possible insight into the critical phenomena, some exact toy models in three 
dimensions were also considered jH], [H]. 

Critical collapse within the scalar-tensor theories has attracted some interest, too 
JU],JT]. There were also some analytical works on the critical collapse in Brans- 
Dicke theory and the theory with conformal coupling J3],521,C3- Using conformally 
transformed Roberts's solutions the authors were able to find analytically the black 
hole mass power law and to extract the critical exponent which turns out to be 7 = 
1/2 ± l/2y2u; + 3 and 7 = 0.21 for the Brans-Dicke and the theory with conformal 
coupling, respectively. 

However, there are serious omissions in the works treating analytically critical col- 
lapse in scalar-tensor theories. First, only two theories were considered. Second, the 
authors have considered simplified models where the matter fields are absent (dilaton- 
vacuum solutions) i.e. the collapsing matter has been taken to be the dilaton field. The 
aim of the present work is to fill in this gap studying analytically the critical collapse 
of a massless real scalar field within the scalar-tensor theories of gravity. Our study 
is based on exact solutions which generalize Roberts's solution in the presence of a 
dilaton field. 

2 Exact self-similar solutions 

The general form of the extended gravitational action in scalar-tensor theories is 

S = f d A xJ^~g (F($)R - Z^)~g^d^d^ 

lo7rG* J v v 

-2U{^))+S m [^ m -g^\. (2) 

Here, is the bare gravitational constant, R is the Ricci scalar curvature with 
respect to the space-time metric g^ u . The dynamics of the scalar field $ depends on 

1 It turns out that the original Roberts's solution as presented in [3] solves the field equations 
only for one nontrivial case and represents just a class of measure zero among the solutions to the 
spherically symmetric homothetic Einstein-scalar field equations given in 0j, |S] and 0- For detail 
discussion see [7| where the Roberts's solution is generalized. I would like to thank L. Burko for 
pointing me this out. 
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the functions F(<&), Z(<&) and U(<&). In order for the gravitons to carry positive energy 
the function F(<&) must be positive. The nonnegativity of the energy of the dilaton 
field requires that 2F($)Z($) + 3[dF($)/d&] 2 > 0. The action of matter depends on 
the material fields \& m and the space-time metric g^ v . 

However, it is much more convenient from a mathematical point of view to analyze 
the scalar-tensor theories with respect to the conformally related Einstein frame given 
by the metric: 

9flv = (3) 
Further, let us introduce the scalar field <p (the so called dilaton) via the equation 

/^\ 2 _3/ rfln(F($)) \ 2 Z($) 

\d$ / 4 y ) 2F($) 1 ' 

and define 

A(<p) = F- 1 ' 2 ^) , 2V(<p) = f/($)F- 2 ($). (5) 
In the Einstein frame action (J2J) takes the form 

S = — ^ J d 4 x^9 (R - 2g^d^d uV - 4V(<p)) 

+S m [V m ,A 2 {<p)g^] ( 6 ) 

where R is the Ricci scalar curvature with respect to the Einstein metric g^ v . 
In what follows we shall consider the case U($>) = and the matter source will be 
a massless real scalar field a. 

The starting point in our considerations is the Roberts's solution which is written 

as 



ds 2 = —dudv + r 2 (u,v)dQ 2 , (7) 

: i ln ^-(i+pK 

2 y u — (1 — p)v 



, 1 , (u - (l+p)v\ 
a = ±-ln ) p , (8) 



where 

r 2 (u, i>) = -[u — (1 — p)v][u — (1 + p)v] 

and p > is an arbitrary constant. 

Dilaton-vacuum solutions can be easily obtained via the conformal transformation 



F \$ dv (u,v)) = A 2 (<p)\ v=t7oiUiV) , (9) 
ds 2 dv = F- l ($ dv (u,v))ds 2 . (10) 

This simple procedure, however, is not applicable when a matter source is present. 
In order to construct exact scalar-tensor solutions describing a collapsing scalar field 
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we shall employ the solution generating methods developed in jHj . Roberts's solution 
is taken seed solution. 

The general form of the scalar-tensor solutions is given by 



l Mu,v)} = F~ 1 [^(a (u,v))] — A 2 [<f(a (u,v))], (11) 
ds 2 = F- 1 [$(<r )]tZs 2 = -F~ 1 [^(a )]dudv + r 2 (u,v)dfl 2 , (12) 
a(u,v) = a((T {u,v)), (13) 



where 



F- l [<&{a G )}r 2 {u,v) (14) 

and the explicit form of the functions F _1 ($(cr )) and a(a ) depends on the partic- 
ular scalar-tensor theory. Below we consider several examples of scalar-tensor theories 
which qualitatively cover the general case and have solutions representable in a closed 
analytical form. 

The solutions above are continuously self-similar since they admit the homothetic 
Killing vector 

d d 

i = u Tu +v Tv- (15) 

The Ricci scalar curvature R is given by 

R = F($(a )) [ 1 + ^ In [F($(a ))] - 1 (J- In [F($(<r ))] ] ) R (16) 
where 

p _ P 2uv n 7 ) 

r ~2^m (17) 

is the Ricci scalar curvature for the Roberts's solution. 

For all solutions we present here the conformal factor F _1 ($) satisfies F _1 ($) > 
A 2 > where A is a constant. Therefore, the conformal factor is everywhere regular 
except at r(u, v ) = where the only singularity 2 of the Ricci scalar curvature is located. 
Thus we can conclude that the causal structure of the scalar-tensor solutions is the same 
as the Roberts's solution. The singularity r{u, v) = is given by u — (1 — p)v and 
u = (1 +p)v for v > and v < 0, respectively. For < p < 1 the singularity is timelike 
while it is null for p — 1. When p > 1 the singularity is spacelike in the region v > 
and timelike for v < 0. 

In all cases we consider, an apparent horizon exists only for p > 1 and it surrounds 
the spacelike singularity while the timelike singularity is naked. In order to get rid of 
the naked singularity we shall proceed as in the case of the Roberts's solution. We can 
replace the region v < by a flat spacetime smoothly matched to the region v > 



2 For all scalar-tensor solutions presented here it can be checked that R could become singular only 
for R. 
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across the v = line where the energy flux vanishes in both the Jardan and Einstein 
frame. 

The fact that the spacelike singularity is hidden behind and apparent horizon signals 
the formation of a black hole. That is why we shall consider bellow the case p > 1 in 
the limit p — > 1. 

With regard to the case < p < 1, we can, just as for the Roberts's solutoin, remove 
the timelike singularity by replacing the region u > and v > by a flat spacetime. 
For < p < 1 the initial conditions are not strong enough to form a black hole. 

The above considerations show that the parameter p is a control parameter and its 
critical value is p cr = 1. Thus the presented solutions can be divided into three classes: 
subcritical (0 < p < 1) , critical (p = 1) and supercritical (p > 1). 



2.1 Brans-Dicke theory 

Brans-Dicke theory is described by the functions = $ and Z(Q) = u/<& where u 

is a constant. Here we shall consider the cases with uj > 0. 

Solution generating transformations pi] give the following scalar-tensor solution: 



•-'<"■•> - ^{M;)- (18) 

where < A 2 < oo. 

The apparent horizon exist only for p > 1 and is located at 

-~ 2(1 ->±v. (20) 



Here and below, in order to present the results in compact and clear form the loca- 
tion of the apparent horizon is determined for (p— 1) << 1. The effective gravitational 
mass inside the apparent horizon is 



1~ ^.J 1 ^ 



1/2 



Mah = - o r AH = ^|A| I — ^ (p - l)V3-VV3+2^. (21) 

As well know the self-similar spacetimes are not asymptotically flat and the black 
hole mass diverges when v — > oo. The standard way to get rid of this problem is to 
truncate the spacetime at some v * > and to add an asymptotically flat region. We 
shall accept from now on that this has been done. 

What is important is that the black hole mass obeys the power law. Hence the 
analytically calculated critical exponent is found to be 7 = | — 2 ^ +2uJ ■ In contrast 
to the dilaton- vacuum case where 7^ = 1 ± 2 ^ 3 1 +2faJ P"T] - [T3] , in the presence of a 
massless scalar field as a matter source, only the minus sign is allowed. 
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2.2 Theory with conformal coupling 

This theory is described by the functions = 1 — |$ 2 and Z($) = 1. The explicit 

functions in the case under consideration are: 



F-\^(u,v)) = ! + (!- A 2 )sinh 2 ^ 



2 / <Tq[U,V) 



(22) 



a(u, v) = VS tanh^ 1 f A tanh j , (23) 

where the parameter A is restricted by < A 2 < 1. 
The apparent horizon exists for p > 1 and is located at 

The mass inside the apparent horizon is 

M AH - ^ (l - A 2 ) V2 V2 (p - l) 1 / 2 - 1 ^,,. (25) 

The critical exponent is 7 = ~ — and coincides with the previously calculated 
critical exponent 7^ ~ 0.21 for the dilaton- vacuum solutions [T2] . 

2.3 F($) = $ and Z($) = ($ 2 - 3$ + 3)/2$($ - 1) theory 
The explicit form the functions $(o"o) and <j(<j q ) is given by 



A 2 

<&~ x (u,v) = — -T7 . y-- — -. — , (26) 

v ; A 2 + (l-A 2 )sm 2 (Aa (u,w)) v 1 

1 + A 2 1 — A 2 
o{u,v) = 2A °o(u,v) — -sin(2Acr (w,f)). (27) 

The range of the parameter A is < A 2 < 1. 
The apparent horizon is located at 



where 



2(1 -p) , \ 

U ^ l±D^p) V <»> 



D xi p) = 1 A(l-A 2 )sin(Aln( P -l)) 
W 2A 2 + (l-A 2 )sin 2 (§Aln(p-l)' 



and D\(p) < 1. 

The mass inside the apparent horizon is given by 
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1/2 

. , n , w (P " l) l7 V (30) 

A 2 + (1 - A 2 ) sin 2 (|Aln(p - 1)) V 1 ± ^P' / 



2 1/2 |A| flfDA(p) 

'AH" = 7 



We see something interesting. The power law M = const(p — p^) 1 observed within 
the general relativity is slightly modified: M = j \p — p cr )(p — Per) 1 ' ■ It is also seen that 
the black hole mass exhibits " damping oscillations" in the control parameter p when 

P-l- 

If we define the critical exponent as 7 = lim^! ln(M)/ ln(p— 1) we obtain 7 = 1/2 
for the case under consideration. 

It is also interesting to find the black hole mass law for the dilaton-vacuum collapse. 
Using solution generating formulas (JHJ) one can show that 

1 (p-l) 1 / 2 

Hence we find 7^ = limp^x ln(M)/ ln(p — 1) = 1/2. 

We have checked that the black hole mass for the near critical collapse in the 
Barker's theory (F($) = $ and Z(<&) = (4 - 3<I>)/2<I>($ - 1)) exhibits "damping 
oscillations" in the critical parameter, too. Since the Barker's case is similar to the 
case under consideration we do not present detail calculations. 

2.4 F($) = $ and Z($) = (1 - 3a 2 $)/2a 2 $ 2 theory 

Here the parameter a satisfies a > 0. 

In the case under consideration we have 



$- % (u,v) = (l + aa (u,v)) 2 +\ 2 , (32) 

a(u, v ) = — arcsin [ I . (33) 

\^(l + aa (u,v)y + X 2 



(j 



Here, A runs < A 2 < 00. 

The location of the apparent horizon is given by 

u = 2(l-p)v. (34) 
The mass inside the apparent horizon is 

Mah = -^a \ ln(p - 1) | (p - 1) V V (35) 

Hence, for the critical exponent we have 7 = lim p ^i ln(M)/ ln(p — 1) = 1/2. 
Respectively, for the dilaton-vacuum collapse we obtain 
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M% = | ln(p - 1) I (p - l) 1 ^ (36) 
and therefore 7^ = 7 = 1/2. 

3 Conclusion 

In the present work we have generalized the general relativistic Roberts-Brady-Oshiro- 
Nakamura-Tomimatsu (RBONT) model in the case of scalar-tensor theories of gravity. 
We have presented new exact self-similar solutions describing a collapse of a massless 
scalar field in scalar-tensor theories. The solutions exhibit a type of critical behavior 
discussed by Choptuik. Three possible evolutions are distinguished: subcritical, critical 
and supercritical. For supercritical evolution the black hole mass low has been found 
and the critical exponent has been extract for several scalar- tensor theories which 
cover qualitatively the general case. It is interesting to note that for some scalar- 
tensor theories within the analytical models we consider, the black hole mass law in 
the near critical collapse is of the form M = f{p — Pcr){p — Per) 1 where f(p — p cr ) is 
a function depending on the particular scalar-tensor theory. One sees that this law is 
slightly different compared to the general relativistic version M = const(p — p cr ) 7 . It 
is also interesting to note that, in some cases, the black hole mass exhibits "damping 
oscillations" in the limit p — > p cr . 

As the RBONT model our models can not be considered as completely realistic 
but they are satisfactory enough in explaining analytically some basic features of the 
critical collapse. 
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